ADA  0  858  64 


TR  79065 


UNLIMITED 


S'  OTIC 

ELF-T - 

FEB  2  8 1980  J 


ROYAL  AIRCRAFT  ESTABLISHMENT  - 


Technical  Report  79065 

June  1979 

ON  THE  USE  IN  STABILITY  ANALYSIS  OF  THE 
RELATIONSHIPS  BETWEEN  THE  COEFFICIENTS 
OF  A  QUARTIC  EQUATION  AND  THOSE  OF 
A  QUADRATIC  FACTOR 


H.H.B.M.  Thomas 


Procurement  Executive,  Ministry  of  Defence 
Farnborough,  Hants 


Thia  document  hew  been  appro? 
for  public  r*’’  t-*  and  tale;  Us 
distribution  is  unlimited. 


2  25  113 


UDC  533.6.013.4  :  512.393  :  517.912.2 


(Bjj 


i+l 


FI 


ROYAL  AIRCRAFT  ESTABLISHMENT 


© 


Technical 


^pmrc-j\ 


9065 


Received  for  printing  7  June  1979 


G 


JON  THE  USE  IN  JTABILITY  ANALYSIS  OF  THE  RELATIONSHIPS  BETWEEN  THE  jgOEFFICIENTS 
OF  A  QUARTIC  ^EgUAT ION  AND  THOSE  OF  A  QUADRATIC  JACTOR  -  / 

by 


"I 


The  stability  of  the  perturbed  longitudinal  motion  of  an  unaugmented  air¬ 
craft  about  a  given  trimned  flight  condition  is  determined  by  the  well-known 
stability  quartic.  It  is  shown  that  two  basic  relationships  exist  between  the 
coefficients  of  the  quartic  and  those  of  a  quadratic  factor.  These  relationships 
are  linear  in  the  coefficients  of  the  quartic  and  the  latter  are  in  turn  linear 
functions  of  the  more  important  stability  parameters.  It  is  possible,  therefore, 
to  use  the  relationships,  already  mentioned,  to  derive  a  further  equation 
connecting  the  two  coefficients  of  the  quadratic  factor  as  a  single  stability 
parameter  is  varied,  all  others  being  held  constant. 


Examination  of  the  position  of  the  locus  traced  out  by  the  point  defined  by 
the  coefficients  of  the  quadratic  in  relation  to  the  axes  and  the  parabolic 
boundary  (separating  real  and  complex  roots)  enables  the  trends  in  the  nature  of 
and  the  degree  of  stability  of  the  mode  or  modes  represented  by  the  quadratic 
factor  to  be  determined.  Expressions  are  given  which  enable  the  value  of  the 
parameter  being  varied,  and  associated  with  a  given  point  of  the  locus,  to  be 
calculated. 


The  use  of  this  (as  far  as  is  known,  novel)  method  of  analysing  aircraft 
stability  is  illustrated  by  a  number  of  examples. 
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INTRODUCTION 


1  2 

It  is  about  a  century  ago  since  Routh  *  formulated  his  well-known 

necessary  and  sufficient  conditions  for  the  stability  of  a  system  as  a  set  of 

inequalities  to  be  satisfied  by  the  coefficients  of  the  characteristic  or 

stability  equation  of  the  system.  Since  the  earliest  days  of  aviation  use  has 

3 

been  made  of  Routh' s  test  functions  (or  their  equivalents  as  given  by  Hurwitz 

A  5 

and  independently  by  Frazer  *  )  to  ascertain  the  stability  or  otherwise  of  a 
given  equilibrium  motion  of  aeroplanes.  However  Routh' s  discriminant  is  a  com¬ 
plicated  expression  even  for  the  quartic  stability  equation,  which  results  from 
the  analysis  of  the  simple  case  of  the  linearized  motion  of  an  aeroplane  with 
motivators  (control  surfaces)  fixed.  It  is  not  surprising,  therefore,  that 
attempts  have  been  made  to  draw  conclusions  about  trends  in  an  aeroplane's 
stability  on  the  basis  of  a  less  complete  knowledge  of  the  coefficients  of  the 
characteristic  equation.  Gates, through  a  series  of  papers  culminating  in  the 
well-known  generalized  theory  of  stability  and  control  analysis  (written 

jointly  with  Miss  H.M.  Lyon)  investigated  the  special  significance  of  the 

2  .... 

constant  term  and  the  coefficient  of  X  of  the  longitudinal  stability  quartic, 
their  relationship  to  the  static  and  manoeuvre  margins  and  the  implications  of 
changes  in  sign  of  these  margins  for  the  stability  of  the  longitudinal  modes  of 
motion  of  an  aeroplane. 

Algebraic  criteria  such  as  those  given  by  Routh  and  Hurwitz  often  lead  to 
conditions  the  interpretation  of  which  is  difficult.  To  overcome  this 
difficulty  and  to  extract  the  utmost  generality  from  specific  calculations  of 
the  roots  of  the  stability,  the  alternative  approach  of  displaying  in  graphical 
form  the  dependence  of  the  stability  on  the  values  of  system  parameters  was 
adopted.  In  the  case  of  a  single  system  parameter,  the  variation  of  the  roots 
of  the  stability  equation  can  be  presented  in  the  form  of  a  one-parameter 
stability  diagram  or  a  root-locus  diagram. 

No  doubt  one  of  the  factors  contributing  to  the  pursuit  of  simplified 
criteria  and  studies  of  general  trends  was  the  tedium  then  attending  the  deter¬ 
mination  of  the  roots  of  the  characteristic  equation.  An  alternative  approach 
was  to  seek  approximate  solutions^  for  the  roots  of  the  characteristic  equation 
Nowadays,  the  relative  ease  with  which  the  roots  can  be  extracted  using  even 
the  small,  readily  available,  desk  computers  could  be  taken  to  indicate  that 
studies  of  general  trends  and  simplified  criteria  become  less  important.  The 
validity  of  this  argument  is  questionable  and  adoption  of  that  viewpoint  could 


4 


result  in  decreased  insight  into  the  nature  of  the  basic  aeroplane  stability. 
Previously  acquired  knowledge  could  become  forgotten  or  at  best  ill-remembered. 
Force  is  lent  to  this  counter  argument  by  the  recent  appearance  of  papers  which 
sought  to  discuss  the  question  of  when  static  stability  is  a  sufficient,  as  well 

g 

as  a  necessary,  condition  for  dynamic  stability.  The  first  of  these  by  Sachs 

investigated  the  matter  using  two-parameter  stability  boundaries  and  root-locus 

9 

plots  for  specific  aeroplanes.  More  recently  Babister  suggested  that  the  use 
of  the  approximate  factorization  of  the  longitudinal  stability  quartic,  due  to 
Bairstow,  afforded  greater  insight  into  the  matter.  He  was  careful,  however,  to 
remark  that  the  conclusions  can  only  be  relied  upon  when  the  degree  of  approxi¬ 
mation  is  good. 

The  quartic  equation  can  be  written  as, 

K,A4  +  K.A3  +  K-A2  4  K,  A  +  Krt  =  0 

5  L  I  U 

or 

K4(A2  +  gA  4  h)(A2  4  GA  4  H)  =  0  , 

where  the  two  factors  always  exist  and  have  real  coefficients.  Previous 
experience  of  the  relationships  between  the  coefficients  of  the  original  quartic 
(K^,  K^,  K2>  Kj  and  KQ)  and  those  of  the  two  quadratic  factors  (g,  h,  G  and  H) 
suggested  to  the  present  author  that  g(G)  is  related  to  h(H)  when  one  para¬ 
meter  of  the  system  is  varied  all  others  being  held  constant.  The  present 
analysis  derives  such  relationships  and  demonstrates  that  they  are  relatively 
simple  in  form.  The  existence  of  these  relationships  enables  diagrams  to  be  con¬ 
structed  from  which  conclusions  may  be  drawn  concerning  the  variation  in  the 
nature  and  stability  of  the  individual  modes  as  a  particular  system  parameter 
changes. 

Although  differing  from  the  one-parameter  stability  diagram  and  the  exact 
form  of  the  root-locus  plot  in  that  it  does  not  require  the  roots  of  the  quartic 
to  be  determined,  the  present  method  provides  all  the  data  necessary  to  construct 
either  of  the  other  two  diagrams,  but  to  do  so  requires  the  solution  of  the 
two  quadratics. 

2  RELATIONSHIP  OF  THE  COEFFICIENT  OF  THE  QUADRATIC  FACTORS  TO  THOSE  OF  THE 
QUARTIC 

In  the  main  text  the  notation  of  R&M  3562  *  is  adopted,  but  for  ease  of 
comparison  with  Babister' s  analysis^  and  the  results  of  other  references**  *2’ ** 
the  analysis  is  repeated  in  Appendix  A  in  terms  of  the  Neumark-type  notation  . 
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The  stability  of  the  small  perturbation  longitudinal  motion  of  an  aircraft, 
with  controls  fixed,  from  a  trimmed  rectilinear  flight  condition  in  a  constant 
density  atmosphere  is  governed  by  the  quartic'® 

K4X4  +  K3X3  +  K2X2  +  K,X  +  K0  =  0  (1) 


where  K. 


and  QmA 


Qm4  » 

+  +  ®wHn3  » 

Qm2  +  ®q0m3  +  ®w^m2  +  mw^m3  +  mu^m3  * 

IDqQm2  +  ®6r^nj  +  ®wHn2  +  ®u^®2  * 

®wHnj  +  ®u^m]  » 

1  +  z*  *  1  , 

Xy  +  zw  +  -  zux^)  *  xu  +  zw  , 

XUZW  -  zuxw  » 

Ue  -  Zq  Ue  , 

ueXu  +  wezu  +  (zuXq  -  xuzq)  -  g2  88  uexu  +  wezu  -  g2  , 

8lzu  ”  82*u  * 

-  xq  -  we  +  (ZqX^  -  z^Xq)  -  uex*  -  wez*  %  -  we  , 

(xwZq  -  zwxq)  -  g,(l  +  z*)  -  g2x*  -  uexy  -  wezu  *  -  g,  -  uexw  -  w£zu 

-  81  zw  +  82^,  • 


The  characteristic  equation  as  given  applies  to  any  body  system  of  axes,  but  it 
is  possible  to  simplify  the  expressions,  without  any  essential  loss  of 
generality,  by  adopting  aerodynamic  body  axes  (so-called  stability  axes)  and 
restricting  our  attention  to  a  level  flight  equilibrium  condition.  In  this  case 
the  following  relationships  apply, 

ug  -  1  ,  wg  -  0  ,  gj  -  g  and  g2  -  0 

It  is  convenient  to  conduct  the  analysis  in  terms  of  dynamic  normalized  equiva¬ 
lents  of  the  physical  quantities.  The  previous  expressions  for  the  W's  and  the 
U's  reduce  to 

W  -  1  ;  U  -  0  ; 

m3  m3 

m2  u  m2  °  w 

W  ■  gz  ;  U  -  -  gz 

m,  ®  u  *  m.  °  w 


6 


so  Chat  Che  coefficients  of  cne  quart ic  may  be  written  as 


Q_  +  m  +  m* 
m~  q  w 


1  m3  q  w 

K_  *  Q_  +  mQ_  +  m*W  +  m 

2  m„  qi,  w  m„  w 


K  *  m  Q_  m*W  +  m  W  +  m  U 
1  q  w  m.  w  m_  u  m„ 


K-  «  m  W  +  m  U 
0  w  m j  u  ni| 

As  remarked  in  the  Introduction  the  quartic  can  always  be  factorized  into  two 
quadratic  factors.  Thus, 

A4  +  K3A3  +  K2A2  +  KjA  +  Kq  -  (A2  +  gA  +  h) (A2  +  GA  +  H)  *  0  (2) 

An  advantage  of  the  two  quadratic  factors  lies  in  the  simpler  statement  for  the 
necessary  and  sufficient  conditions  for  stability  in  terms  of  g,  h,  G  and  H 
rather  than  in  terms  of  K^,  K2,  Kj  and  .  For  the  former  the  conditions 
are: 

g>0,  h  >  0  ;  G  >  0  ,  H  >  0 

whilst  in  terms  of  the  coefficients  of  the  original  quartic  the  necessary  and 
sufficient  conditions  for  stability  are 

K3  >  0  ,  Kj  >  0  ,  KQ  >  0  and  R  >  0 

2  2 

where  R  (the  Routh  discriminant)  is  -  Kj  -  K-jKq  • 

An  essential  step  in  the  analysis  in  terms  of  g,  h  or  G,  H  is  the 
relationship  between  these  and  the  coefficients  K3,  K2>  Kj  and  Kq  . 

Expansion  of  the  product  of  the  two  quadratic  factor  yields  the  following 


G  +  g 


H  +  h  +  gG  ■ 
Hg  +  Gh  - 


i 
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Since  subsequently  we  shall  regard  and  as  parameters  which  may  be 

varied,  whilst  other  quantities  involved  in  K^,  K2,  Kj  and  Kq  are  held 
constant,  the  coefficients  are  rewritten  as 


K.  =  k  +  m 
2  2  w 


K.  =k,+mW  +mU 
1  1  w  m2  u  m2 


K_  =  m  W  +  m  U 
0  w  mj  u  ni| 


Manipulation  of  the  above  relationships  between  the  two  sets  of  coefficients 
yields  the  following  equations, 


K2g  -  Kj  +  K3h  =  2gh  +  g  (K3  -  g) 


K2h  -  Kq  =  gh(K3  -  g)  +  ti 


An  identical  pair  of  equations  relate  G  and  H  by  virtue  of  the  synmetry  of 
the  basic  equations. 

9 

To  obtain  the  Bair stow  approximations  used  by  Babister  the  assumption  is 
made  that  g  and  h  ,  which  may  be  identified  with  the  phugoid  motion,  are 
small  so  that  to  this  degree  of  approximation. 


K2h  “  K0  *  0  °r  h  =  K^ 


K2g  -  K,  +  K3h  =  0  or  g 


<K2K|  -  k3k0> 


Furthermore  since  g  and  h  are  small. 


G  -  K3  -  g  «  K3 


H  =  K2  -  h  -  gG  K„ 


These  various  approximations  may  be  improved  upon  by  an  iterative  solution.  How¬ 
ever  our  objective  is  not  to  obtain  approximate  values  of  g  and  h  or  the 
roots  associated  with  this  factor,  but  to  examine  the  manner  in  which  g  and  h 


are  related  when  a  single  parameter  (for  example  mu>  is  varied.  We  accordingly 
return  to  equations  (4). 

3  RELATIONSHIP  OF  g(G)  TO  h(H)  WHEN  mu  IS  VARIED 

Introduction  of  the  expressions  for  Kj  and  Kq  into  equations  (4) 
yields  the  following  pair  of  equations, 


U  m  =  K.g  +  K_h  -  k.  -  Wm  (K  -  k_)  -  2gh  -  g^(K-  -  g) 
m2  u  z  j  i  m2  z  z  j 

U  m  -  K,h  -  W  (K.  -  k,)  -  gh(K  -  g)  -  h2  . 
ni|  u  z  nij  z  z  j 


These  two  equations  imply  that  as  m^  is  varied  the  associated  values  of  g 
and  h  satisfy  the  relationship, 

U  k,  +  W  (K„  -  k0)  -  K_h  -  K.g  +  2gh  +  gZ (K_  -  g) 


m,  Wm  (K2  -  k2)  -  K2h  +  gh(K3  -  g)  +  h 

This  exact  relationship  between  g  and  h  can  be  rearranged  as  a  quadratic  for 
h  in  terms  g  ,  thus, 

a2h2  +  ajh  +  aQ  =  0  (7) 


where  a2  =  k 


a,  =  -  kg^  +  (K3k  -  2)g  +  (K3  -  K2k) 

«0  =  g3  “  k382  +  K28  +  (kWml  "  Wm2)(K2  "  V  "  kl  • 

Approximations  to  equation  (7)  may  be  formed. 

(1)  Treat  both  g  and  h  as  small  (appropriate  when  isolating  the  phugoid 
factor) 

Equation  (7)  may  be  then  approximated  by  the  linear  relationship 


“lh  +  aQ  "  0 


where  ctj  -*■  (K3k  -  2)g  +  (K3  -  K2k) 

and  aQ  -+  K2g  +  (k%j  -  Wm2)(K2  -  k2)  -  kj  . 

9 

This  is  the  same  result  as  obtained  by  Babister  by  use  of  the  Bairstow  approxi¬ 
mation  to  the  roots  of  the  quartic. 
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(2)  Retain  terms  up  to  and  including  the  product  and  the  squares  of  g  and  h 

The  equation  has  the  same  form  as  (7)  with  01q  coefficient  simplified  by 
the  omission  of  the  term  g3  ,  an  often  insignificant  change. 

To  each  pair  of  values  of  g  and  h  that  satisfy  equation  (7)  there 
corresponds  some  value  of  the  parameter  mu  and  this  is  given  by  either  of  the 
two  relationships  of  equation  (5). 

As  previously  remarked  a  corresponding  set  of  equations  exists  between 
the  other  two  coefficients  G  and  H  .  In  a  complete  analysis  it  is  useful, 
as  we  shall  see  later,  to  plot  the  curve  of  H  as  a  function  of  G  as  well  as 
of  h  as  a  function  of  g  ,  because  in  such  a  case  extreme  values  of  the  design 
parameter  need  to  be  embraced.  Where,  on  the  contrary,  interest  centres  around 
the  effect  of  modest  variations  from  a  given  datum  value  it  is  more  reasonable 
to  obtain  either  approximate  or  exact  values  of  g  and  h  for  the  basic 
aircraft. 

To  illustrate  the  interpretation  of  the  plot  of  h  against  g  this  has 

been  displayed  schematically  in  Fig  1.  As  previously  noted  stability  is  ensured 

if  g  >  0  and  h  >  0  .  This  positive  quadrant  is  subdivided  by  the  curve 
2 

g  *  4h  into  regions  in  which  the  factor  corresponds  to  a  convergent  oscilla- 

2 

tion  and  two  subsidences  respectively.  The  axes  and  the  curve  g  =  4h  sub¬ 
divide  the  unstable  regime  into  different  forms  of  instability  as  indicated  on 
the  diagram. 

3. 1  Value  of  g  when  h  =  0 

Of  special  interest  is  the  value  of  g  when  h  is  zero.  The  determina¬ 
tion  of  this  value  of  g  requires  the  solution  of  a  cubic  equation  =  0  , 
that  is 

g3  -  K3g2  +  K2g  -  k,  +  (K3  -  k2)(kWm]  -  Wm2)  =  0  . 

The  different  roots  corresponds  to  the  'phugoid'  factor,  which  we  shall 

generally  identify  with  g  and  h  ,  and  to  the  'short  period'  (or  fast)  factor 

3  2 

(G  and  H) ,  In  the  case  of  'phugoid'  factor  we  may  drop  the  term  in  g  and  g 
for  the  numerically  smallest  root.  It  is,  therefore,  of  interest  to  examine 
mor’;  closely  the  constant  term  in  the  cubic  for  g  ,  which  we  shall  denote  by 
K  .  Now 


I 
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K  =  (K3  -  k2) (kWmi  -  Wm2)  -  k, 

=  ®  ,{(g  +  x  )  —  -  x  >  -  m  (x  z  -  x  z  )  -  gz  m. 
wj  wzw  uf  q  uw  w  u'  6uw 

If,  for  simplicity,  we  neglect  the  effects  of  compressibility  and  elastic  distor¬ 
tion  the  dynamic  normalized  quantities  involved  in  the  above  expression  can  be 
rewritten  in  the  following  forms, 

8  =  CL. 


x  =  -  X  =  2 (C_.  +  C.  J 

u  u  D  AS'e  ’ 


where  C 


AS 


.  3T 

1  av 

~  ~  3Q  »  so  that  CAS  =  0  for  jet  propelled  aircraft, 
-2  =  2Ct 


w 


w 


g  +  x 


'  K  •  (t?  ■  cl)6 
-  -  (t£ ♦  %\  *  (^)e 

(?\\  r%\ 

\  4  ‘  1 4W4  ' 


and  so 


The  first  term  in  K  is  thus  approximately. 


'  2mw(CD  *  CAS  "  CL  3C  ) 


which  demonstrates  the  importance  of  whether  the  equilibrium  speed  is  above  or 
below  the  minimum  drag  speed  (c/  Refs  7,  9,  11)  in  determining  the  sign  of  g 
when  h  =  0  . 

If  the  basic  value  of  g  is  positive  there  are  three  cases  according  as. 


(i)  (g), 


»  0  ; 


S90fc 
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(iii)  (g) 


basic 


(8) 


h=0 


as  illustrated  in  Fig  1.  Case  (iii;  arises  when  K  is  positive,  that  is,  when 
(CD  +  CAg  -  Cl(8Cd/8Cl))  has  a  sufficiently  large  negative  value.  When  this  is  so 
there  is  a  range  of  positive  static  margins  for  which  undamped  phugoid  oscilla¬ 
tions  can  occur  (of  Appendix  A  and  examples).  An  examination  of  the  remaining 

terms  in  K  indicates  that  small  values  of  m  and  m*  make  for  more  positive 

q  w  r 

values  of  K  and  the  same  consequences,  which  is  in  agreement  with  the  conclu- 

.  9 

sions  drawn  by  Babister  . 

4  THE  g,  h  RELATIONSHIP  FOR  VARYING  m 

w 


Another  important  and  basic  parameter  in  the  longitudinal  stability 
quartic  is  m^  .  The  relationship  of  g  to  h  (or  G  to  H)  as  is 

varied  can  readily  be  found  by  introduction  of  the  expressions  for  K2>  Kj  and 
Kg  into  equations  (4)  and  rearranging  the  equations  obtained  in  the  form. 


mw(g  -  Wm2) 
mw0i  "  Wm)) 


kj  +  Um2mu  -  k2g  -  K^h  +  2gh  +  g2(K3  -  g) 
*  k2h  +  gh(K3  -  g)  +h2 


(8) 


Elimination  of  m^  yields  the  desired  relationship,  that  is 


8  -  Wm2 

h  -  Hn, 


kl  +  Ura2mu 


Un>lmu 


-  k2g  -  K3h  +  2gh  +  g2(K3  -  g) 

-  k2h  +  gh(K3  -  g)  +  h2 


(9) 


Equation  (9)  is  the  most  general  form  of  the  relationship.  It  is  instructive 
to  introduce  into  the  above  equation  the  simplified  forms  of  Umj ,  Um2,  Wm]  and 
Wm2  which  pertain  to  a  level  flight  condition  of  equilibrium.  We  then  obtain 
the  equation, 


8  ~  Xu 

h  '  ZuCLe 


-  mu(CLe  +  xw)  -  k2g  -  K3h  +  2gh  +  g  (K3  -  g) 
h2  +  gh(K3  -  g)  -  k2h  -  zwmuCLe 


(in) 


This  may  be  rewritten  as  a  quadratic  in  h  , 


+ 


S,h 


0 


(11) 
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where 


ik2Xu  ‘  kl  +  mu(CLe  +  Xw}  ‘  ZuCLeK3*  +  (2zuCLe  “  XuK3)g  +  Xug2  * 

CLe[lklzu  +  Vumu  “  Zumu(CLe  +  xw^  ”  (k2Zu  +  zwmu)g 

+  Zug2(K3  “  g)]  ' 


From  this  equation  the  values  of  h  may  be  obtained  as  g  is  varied  away  from 
its  basic  value  as  shown  schematically  in  Fig  2.  The  interpretation  of  the 
different  regions  into  which  the  axes  and  the  curve  g  =  4h  divide  the  plane 
follows  the  pattern  of  Fig  1 . 

Various  approximations  to  the  relationship  of  h  to  g  may  be  obtained. 
The  simplest  of  these  corresponds  to  the  retention  of  only  linear  terms  on  the 
right-hand  side  of  equations  (8)  or  in  the  numerator  and  denominator  of  the 
right-hand  side  of  equation  (10).  To  this  approximation  (approximation  I) 


B2  *  K3  • 


0  88  k  x  +  m  (Ct  +x)  -  k,  ■  z  Ci  L  , 

1  2  u  u  ^e  w  1  u  Le  3 

0n  155  CT  jk  z  +  x  z  m  -  m  z  (Ct  +  x  )}  -  (z  m  +  k„z  )g 

0  Le  I*  1  u  uwu  uu^e  w  ’  wu  2  u  ° 


Another  form  of  approximation  is  obtained  if  in  the  expanded  equation  (11) 
2  2  3 

terms  in  gh  ,  g  h  and  g  are  neglected.  With  this  adjustment  (approximation 
II) 

®2  ^3  ~  xu  » 

6.  858  jk„x  -  k  +  m  (Cr  +  x  )  -  z  Cl  K.J  +  (2z  Ct  -  x  K„)g  , 

1  '  2  u  1  u  ^e  w  u  e  3’  u  Le  u  3  ° 

8n  35  Ct  |{k,z  +  x  z  m  -  z  m  (Ct  +  x  )[  -  (k„z  +  z  m  )g  +  z  K_g2| 

0  LeM  1  u  u  u  u  u  u  Le  w  >  2  u  w  u  u  3° 


With  this  degree  of  approximation  the  curve  becomes  a  conic  and  in  Appendix  B 
the  nature  of  conic  approximations  to  the  h,  g  curves  for  both  varying  m^ 
and  varying  mw  are  discussed. 
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It  is  implicit  in  the  above  approximations  (and  those  of  section  4.1)  that 
g  and  h  are  generally  small  and,  in  the  case  of  the  longitudinal  quartic, 
associated  with  the  slow  or  'phugoid*  motion. 

4. 1  Value  of  g  for  which  h  =  0 

The  exact  relationship  between  h  and  g  reduces  to  the  following  cubic 
in  g  when  h  is  set  to  zero. 


3  2 

z  g  -  z  K»g  +  (k_z  +  z  m  )g  -  Ik.z 
u  u  3°  2  u  w  u  °  111 


x  z  m 
u  w  u 


-  z  m  (Ct 
u  u  ^e 


x  )}  = 

w  * 


As  a  first  approximation  to  the  solution  of  this  cubic  we  have. 


g  = 


kz  +  m  lx  z  -  z  (Ct  +x){ 
1  u  ut  u  w  u  ^e  w  * 


k„z  +  z  m 
2  u  w  u 


kl%  -  mu<(CD  *  CAS)(Cl-(,  *  V  - 
Vie  *  ”UCL„ 


where  Ct  *  (3CT /3a) 

L  e 

The  term  in  in  the  numerator  is  approximately  equal  to 


muCLa(CD  +  CAS  “  CL  3Cl) 


Thus  we  see  that  once  more  the  location  of  the  flight  conditions  relative  to  the 
minimum  drag  conditions  plays  an  important  part  in  determining  the  trends  in  g  . 

5  NUMERICAL  EXAMPLE 

To  illustrate  the  use  of  the  preceding  analysis  we  consider  a  basic  air¬ 
craft  having  the  following  characteristics  and  equilibrium  conditions, 


L,  =  272 

iy  ■  '-75 

CLe  -  0.25 

x  =  0.0188 

u 

z 

u 


0.25 
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X 

w 

= 

-0.139 

z 

w 

= 

4.899 

m 

w 

= 

113.774 

m« 

w 

= 

1 .20 

m 

q 

= 

6.542 

and  with  z»  ;  z  and  m  all  zero, 
w  q  u 

To  examine  the  effect  of  varying  m^  on  the  coefficients  g  and  h , 
equation  (7)  is  solved  for  the  values  of  g  indicated  in  Fig  3.  In  this  example 
the  basic  aircraft  possesses  a  stable  phugoid  mode  (convergent  oscillation) . 

As  mu  is  increased  a  region  of  two  subsidences  (very  narrow)  is  quickly  tra¬ 
versed  and  a  region  in  which  the  quadratic  factor  corresponds  to  two  real  roots, 
one  subsidence  and  one  divergence,  is  entered.  Reduction  of  m^  is  associated 
with  a  phugoid  or  oscillatory  motion  at  first  convergent  but  eventually  becoming 
divergent. 

By  use  of  equation  (11)  the  effect  of  varying  mw  may  be  determined. 

The  values  of  h  for  the  values  of  g  indicated  in  Fig  4  yield,  for  the 
example  aircraft,  the  curve  shown  in  that  figure.  As  the  value  of  mw  is 
increased  this  curve  is  traversed  in  the  direction  indicated.  Most  of  the  curve 
remains  within  the  convergent  oscillation  region  and  it  just  enters  the  one  sub¬ 
sidence,  one  divergence  region  near  g  *  0.03  . 

6  VALUES  OF  THE  PARAMETERS  m  AND  m  ASSOCIATED  WITH  POINTS  OF  THE 

_ u _ w _ 

h,  g  CURVE 

To  each  point  of  the  h,  g  curve  there  corresponds  a  value  of  the  design 
parameter  being  varied,  namely,  m^  or  m w  ,  as  the  case  may  be. 

The  values  are  given  by  either  of  the  two  equations  (5)  for  m^  ,  when 
this  is  the  parameter  being  varied  and  by  either  of  the  two  equations  (8)  for 
nv  when  this  is  being  varied.  When  either  (or  both)  of  the  factors 
(A^  +  gA  +  h)  and  (A^  +  GA  +  H)  represents  an  oscillation  no  difficulty  arises 
in  the  computation  of  the  parameter  values.  However,  when  these  two  quadratic 
factors  themselves  factorize  into  two  pairs  of  real  roots  the  matter  is  not 
straightforward.  Suppose  the  factors  are  (A  +  Rj),  (A  +  R2),  (X  +  Tj),  (A  +  v^) 
then  all  six  combinations  of  these  in  pairs  must  occur  somewhere  along  the  g,  h 
curve  or  the  G,  H  curve,  since  the  analysis  cannot  make  the  distinction.  It 
is  for  this  reason  that  the  simultaneous  plotting  of  the  g,  h  and  the  G,  H 
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curve  is  desirable.  The  matter  is  discussed  in  more  detail  in  relation  to  the 
numerical  examples  given  in  Appendix  A. 

7  ANALYSIS  IN  TERMS  OF  THE  PARAMETERS  K  ,  STATIC  MARGIN  AND  H  , 
_ n  _ m 

MANOEUVRE  MARGIN 

The  well-known  criteria  for  stability  of  the  longitudinal  motion  of  an 

aircraft  (static  and  manoeuvre  margins)  are  related  to  m  and  m  as  follows, 

u  w 


K 

n 

and 

H 

m 


v]  =Le<CLa)e  {(CL  +  2  3v)e  mw  "  *(CLa  +  CD)emu} 


— /  '  -  x  (m  +  Ci  m  ) 

1*  i  ^cLa^e  W  a  9  e 


{of  Appendix  A) . 


Thus  the  coefficients  of  the  quartic  may  be  expressed  in  terms  of  K  and  H 

n  m 

in  place  of  m  and  m  .  The  coefficients  affected  are  K„,  K,  and  K_  . 
r  u  w  2  10 

These  may  be  written  in  the  form. 


and 


K2  "  K2mHm+K20 


K,  *  K,  H  +  K,  K  +  K, _ 
1  1mm  In  n  10 


K.  =  Kn  K 
0  On  n 


where  K 

K 

K 

K 

K 

K, 


20 

= 

* 

s°* 

1 

CM 

10 

s 

kl  '  mqiVw  “  Zu(CLe 

‘2m 

= 

Vw;iy  » 

'lm 

- 

(VVlVw  "  zu(CLe 

‘In 

* 

<2u,/iy)CLe(CLe  +  xw) 

On 

- 

(2y./i  )C?  z  , 

1  y  Le  w  ’ 

when  the  approximation 


w 


8CL 

CD 


9C. 


_ L 

9a 


is  made 
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By  introduction  of  the  relationships  appropriate  to  level  flight,  viz. 


the  above  can  be  expressed  in  terms  of  C_,  C.  and  their  derivatives. 

\J  L 


With  the  coefficients  of  the  quartic  expressed  in  terms  of  and 

we  can  establish  the  locus  of  the  points  g,  h  and  G,  H  as  each  parameter  is 
varied  in  turn,  all  other  parameters  being  held  constant.  From  equations  (4) 
we  have, 

'  K,0  '  *20*  '  K3h  *  2<*  *  82<>S  -  8)  ' 
and  (12 

Vm  -  K0„K„  *  h2  *  S»«3  -  8>  -  t20h  • 


Elimination  of  gives  the  locus  for  varying  and  this  is 


(K10  +  KlmV  ‘  (K?n  +  K7m)H„  -  K,h  +  2gh  +  gZ(K,  -  g) 


20  2m  m 


h  +  gh(K3  -  g)  -  h(K20  ♦  K2mHm) 


Elimination  of  yields  the  relationship  satisfied  for  varying  ,  that  is 


K2mS  -  K1m 


(K10  +  KlnV  *  K208  ~  K3h  *  2gh  *  §  (K3  ~  8> 
h2  +  gh(K3  -  g)  -  hK2Q  +  K0nKn 


.  (14) 


By  virtue  of  the  linear  relationship  between  H  and  m  ,  a  fixed  value  of  one 

m  w 

implies  a  fixed  value  of  the  other.  Consequently,  the  curves  for  varying 
(equation  (13))  is  identical  to  equation  (6).  However  a  fixed  value  of  the 
static  margin,  ,  implies  a  linear  relationship  of  my  and  m^  as  is 
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varied.  In  this  case  Kg  maintains  a  constant  value  proportional  to  Kn  . 

This  is  illustrated  by  the  results  shown  in  Figs  18  and  19,  although  it  should 
be  noted  that  these  results  are  in  terms  of  the  old  notation  (see  List  of  symbols 
and  Appendix  A) . 

8  DISCUSSION 

The  analysis  of  a  number  of  different  cases  enables  some  comments  to  be 
made  as  regards  the  stability  of  the  different  longitudinal  modes  of  motion  and 
their  dependence  on  design  features.  These  are  in  agreement  with  those  of 
Refs  6,  7,  8,  9  and  1 1 . 

From  the  equation  of  the  g,  h  locus  for  varying  k  and  varying  u> 
respectively  it  is  possible  to  deduce  the  value  of  g  when  h  ■  0  and  to 
demonstrate  its  strong  dependence  on  the  relation  of  the  equilibrium  speed  to 
the  minimum  drag  speed.  This  is  illustrated  by  Figs  5  and  6. 

Cases  can  exist  for  which  increase  in  the  static  margin  (K^  or  decrease 

in  k)  results  in  an  unstable  phugoid  oscillation,  for  example,  aircraft  3  of 

2 

Fig  5.  It  is  worth  noting  that  the  factor  (X  +  GX  +  H)  represents  the  usual 
short-period  oscillation  of  aircraft  1 ,  2  and  3  and  that  the  values  of  G  and 
H  only  change  slightly  for  the  range  of  g  covered  in  Fig  5.  Each  of  the 
aircraft  1,  la  and  3  of  Figs  5  and  6  is  statically  stable  (E  >  0)  for  points 
within  the  unstable  oscillation  region  (see  also  section  4  of  Appendix  A). 

Variation  of  to  with  all  other  parameters  held  constant  shows  that  air¬ 
craft  I  and  2  (Fig  7)  possess  unstable  phugoids  for  a  certain  range  of  to  in 
spite  of  being  statically  stable  (see  section  5  of  Appendix  A). 

The  cases  illustrated  in  Figs  8  and  9  point  to  the  danger  of  a  positive 
static  margin  in  conjunction  with  a  negative  or  zero  manoeuvre  margin,  which 
can  give  rise  to  a  rapid  divergence  or  a  markedly  unstable  phugoid  oscillation. 

It  should  be  noted  that  for  the  cases  illustrated  in  Figs  8  and  9  when  both  g 
and  h  are  positive,  H  must  be  negative  for  a  negative  static  margin  (K^) 

(or  a  negative  E) .  This  implies  a  divergence  in  pitch.  In  contrast, with  the 
values  of  manoeuvre  margin  indicated  in  Fig  10, a  divergence  is  associated  with 
the  negative  static  margin.  As  pointed  out  in  Ref  6  this  divergence  may  be 
sufficiently  slow  to  be  tolerated  if  the  manoeuvre  margin  is  sufficiently  large 
and  positive  as  suggested  by  the  trend  with  increasing  manoeuvre  margin. 

With  a  fixed  value  of  k  ■  -0.085  ,  variation  of  to  ,  for  the  aircraft 
which  is  the  subject  of  Figs  13  and  14,  demonstrates  stability  of  all  modes  for 


u)  >  -2.663  .  (This  corresponds  to  K  *  0,  H  «•  0.01  .)  Analysis  in  terms 

n  m 

of  manoeuvre  margin  for  the  same  aircraft,  but  with  the  static  margin  held  at 
0.05,  again  demonstrates  the  link  between  unstable  phugoid  oscillation  and  a 
small  manoeuvre  margin  with  this  level  of  static  stability. 

9  CONCLUSIONS 

(1)  It  is  possible  to  relate  the  coefficients  of  the  two  quadratic  factors, 
into  which  it  is  always  possible  to  factorize  a  quartic  stability  equation,  to 
those  of  the  quartic. 

(2)  The  relatively  simple  nature  of  the  relationships  between  the  coefficients 
provides  a  method  of  analysing  the  stability  of  the  longitudinal  motion  of  an 
aircraft  as  a  single  design  parameter  is  varied. 

(3)  A  stability  analysis  of  this  form  is  intermediate  between  the  one-parameter 
diagram  and  a  root-locus  plot.  It  does  not,  however,  involve  solution  of  a 
number  of  quartics.  It  also  possesses  a  compactness  certainly  not  present  in 
the  one-parameter  diagram. 

(4)  Use  of  the  present  method  enables  certain  broad  trends  to  be  established, 
see  section  8. 

(5)  Here  the  application  of  the  method  has  been  to  the  longitudinal  stability 
quartic.  It  is  equally  applicable  to  analysis  of  the  lateral  stability  quartic 
and  extension  to  higher-order  stability  equations  may  be  possible. 
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Appendix  A 

ANALYSIS  OF  THE  QUARTIC  IN  TERMS  OF  THE  NOTATION  OF  R&M  2027 

A.  1  Analysis  for  the  design  parameters  id  and  ic 
9 

Babister  uses  the  notation  of  R&M  2027  and  to  facilitate  comparison  with 
his  analysis  and  the  results  of  a  number  of  past  papers  it  is  useful  to  analyse 
the  quartic  using  this  old  notation. 

The  stability  quartic  for  small  perturbation  longitudinal  motion  about 
steady  level  flight  in  a  uniform  atmosphere  has  the  same  form,  but  the  numerical 
value  of  the  coefficients  and  roots  differ  from  those  of  the  dynamic  normalized 
equation  of  the  main  text.  We  write  it  in  the  form, 


X4  +  BX3  +  CX2  +  DX  +  E 


(A- 1) 


where  (see  Refs  6,  7,  9  or  II) 

B  ■  N  +  v  +  x  » 

C  =  w  +  Nv  +  Qx  +  P, 

D  =  Qu  -  Sk  +  Pv  +  Rx  , 

E  ■  Rio  -  Tic  . 

Corresponding  to  equations  (4)  of  the  main  text,  section  2,  we  have. 


Cg  -  D  +  Bh  -  2gh  +  gz(B  -  g) 


Ch  -  E  =  gh(B  -  g)  +  h 


(A-2) 


If  the  above  expressions  for  C,  D  and  E  are  inserted  in  these  equations 
they  can  be  written  as 


(g  ~  Q)<u  +  Sk  *  Pv  +  Rx  -  (P  +  Nv  +  Qx)g  -  Bh  +  2gh  +  g  (B  -  g) 


(h  -  R )id  +  Tic  *  -  (P  +  Nv  +  Qx)h  +  gh(B  -  g)  +  li 


(A-3) 


Thus,  if  we  require  the  relationships  of  h  and  g  as  ic  is  varied  all  other 
parameters  being  held  constant  it  is  necessary  to  eliminate  ic  .  This  yields 


S  (Pv  +  Rx  +  Q(»)  -  (id  +  P  +  Nv  +  Qx)g  -  Bh  +  2gh  +  g  (B  -  g) 


Rid  -  (id  +  P  +  Nv  +  Qx)h  +  gh(B  -  g)  +  h^ 


.  (A-4) 
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For  varying  w  the  relationship  is  given  by 


8  -  Q  ,  (Pv  ♦  Rx  ~  Sk)  -  (P  +  Nv  ♦  Qx)r  -  Bh  +  2gh  +  g^(B  -  g)  .  0 

h  -  R  ’  *  K 


-  Tk  -  (P  +  Nv  ♦  Qx)h  +  gh(B  -  g)  +  h 


As  remarked  in  the  main  text  the  value  of  <  corresponding  to  a  given  point  of 
the  curve  defined  by  equation  (A-4)  is  given  by  either  of  the  two  equations. 


Sk  -  (Pv  +  RX  +  Qco)  -  (ui  +  P  +  Nv  +  QX)g  -  Bh  +  2gh  +  g^(B  -  g)  (A-6) 


or 


Tic  *  Rw  -  (to  +  P  +  Nv  +  Qx)h  +  gh(B  -  g)  +  h 


(A-7) 


Likewise  for  a  point  on  the  curve  defined  by  equation  (A-5)  the  value  of  a>  is 
given  by  either 

(g  -  Q)ai  -  (Pv  +  Rx  -  Sic)  -  (P  +  Nv  +  Qx)g  -  Bh  +  2gh  +  g2(B  -  g)  (A-8) 
or 

(h  -  R)ai  -  -  Tic  -  (P  +  Nv  +  QX)h  +  gh(B  -  g)  +  h2  .  (A-9) 

It  may  be  noted  that  as 


also  as 


h  -*■  R  ,  w  ■+■  ±  00 

G  -*■  Q  ,  a)  -*•  ±  ® 


Since  equations  (A-8)  and  (A-9)  yield  identical  values  of  to  (this  generally 
provides  an  arithmetical  check)  the  point  (Q,R)  is  a  point  of  the  curve,  which 
approached  from  one  side  corresponds  to  <0  -*•  +»  and  from  the  other  side  to 
u  .  Since  Q  and  R  are  usually  positive  this  point  lies  in  the  first 

quadrant . 

A  discussion  of  the  geometry  of  the  curves  defined  by  equations  (A-4)  and 
(A-5)  is  given  in  Appendix  B. 

A. 2  Analysis  in  terms  of  H  and  K 
_ m _ n 

Instead  of  expressing  the  coefficients  C,  D  and  E  in  terms  of  10  and 

<  we  may  express  them  in  terms  of  the  stability  margins  H  and  K  and  then 

in  n 

examine  the  effect  of  varying  either  of  these  two  parameters  on  the  stability 
characteristics. 
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The  static  and  manoeuvre  margins,  K  and  H  ,  as  usually  defined  are 

*  n  m 


B  1  ,  A  v 

— — c r (w  • 


(A- 10) 


(A- 1 1 ) 


Here  level  flight  equilibrium  conditions  have  been  assumed  and  the  approximation 
J  (CLa  +  Cp)  *  roade.  From  these  relationships  we  have, 

a)  =  JL  (JL)  c,  H  -  JCT  v 
iB  \2H )  La  m  2  La 

and 

*  ■  q  (i)(cL  *  i  ^r)  H»  -  ^  (!)  CLK»  -  (CL  *  i  *-w) “ 


so  that 


C  =  C  H  +  (N  -  £Ct  )v  +  QX  +  P  , 
mm 


D  =  DH  +  DK  +<P  + 
mm  n  n  ] 


E  =  E  K 
n  n 


v  +  Rx  , 


where  C 


t  (n) ’ 

t  (i)fc^  -  (cl  * !  if)  s)  • 

t  (f)  V  • 
t  (*)  c^«  • 


If  the  expressions  for  C,  D  and  E  are  introduced  into  equations  (A-2) ,  we 
have  in  place  of  (A-3)  the  following, 
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(Cg-D)H  -DK  =  jp  +  (C.  +  JV(3Ct/3V))s}v  +  Rx  -  jP  +  Qx  +  (N  -  JCr  )v}g 
m  m  m  n  n  *  L  L  ’  “  ' 

-  Bh  +  2gh  +  g2(B  -  g) 

► 

and 

hCmHm  "  EnKn  =  “  {p  +  Qx  +  (N  -  iCLa)v}h  +  gh(B  -  g)  +  h2  , 


.  (A-12) 

which  yield  for  the  relationship  between  h  and  g  for  varying  only  the 

equation, 

tDmHm  +  RX  +  lp  +  <CL  +  |V(3Cl/3V))s|v]  -  +  P  +  Qx  +  (N  -  iCLa)v}g 

°n  _  _ _ _ -  Bh  +  2gh  +  g2  (B  ~  g) 

En  h2  +  gh(B  -  g)  -  jCmHm  +  P  +  QX  +  (N  -  jCLa>v}h 


.  (A- 13) 

whilst  if  K  is  held  constant  along  with  the  other  parameters  and  H  varied 
n  m 

h  and  g  are  related  according  to  the  equation, 


(C.8  -  V 

C  h 
m 


{dk  +  (P  +  CL  +  }V(3Cl/3V))v  +  Rx}  -  {p  +  Qx  +  (N  -  {cLa)v}g 

-  Bh  +  2gh  +  g2  (B  -  g) 

— i  i  -  -  -  ^  * 

E  K  -  jP  +  Qx  +  (N  -  iCL  )}h  +  gh(B  -  g)  +  h 
n  n  *  ci 


.  (A- 1 4) 

From  the  equations  for  w  and  <  in  terms  of  H  and  K  it  can  be  deduced 

m  n 

that  equation  (A-13)  defines  the  same  curve  as  equation  (A-4) .  In  contrast  vary¬ 
ing  (with  held  constant)  may  not  be  identified  with  varying  o>  ,  since 

it  also  implies  a  varying  value  of  k  .  This  reflects  the  natural  character  of 
u  and  <  (or  K  )  as  parameters  of  the  aircraft  or,  in  other  words,  the  fact 

that  K  is  based  on  an  exact  theory  whilst  H  is  the  result  of  an  approxi- 
n  m 

mate  theory.  The  point  g,  h  corresponds  to  some  value  of  and  and 

these  latter  may  be  deduced  from  either  of  the  two  equations  formed  from 
equations  (A-12)  (of  equations  (A-6)  to  (A-9)).  The  value  of  -*•  ±“  as  the 
point  h  =  0,  g  =  D  /C  is  approached. 
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A. 3  Numerical  examples 

Some  of  the  features  of  the  analysis  are  best  illustrated  by  numerical 
examples.  Again  to  aid  in  comparison  of  the  results  and  conclusions  with  those 
of  earlier  investigations  the  numerical  work  is  in  terms  of  the  old  notation  and 
the  basic  data  drawn  from  Refs  6,  7  and  II. 

The  characteristics  of  the  four  basic  aircraft  considered  initially  are 
listed  below. 


Aircraft 

CL  = 

0.8  , 

0) 

=  25.00 

K  —  0  , 

V  = 

1  ,  x  =  0.4  , 

CLa 

8Cl 

=  1ST  =  4-° 

,  CAS  =  0.01  , 

CD0  = 

0.01  , 

s 

=  0.0625 

°D  ■  \  +  sCL  '  °-05 

Aircraft 

_2 : 

CL  = 

0.2  , 

0) 

=  25.00  , 

k  =  -  21.00  , 

V  = 

1  ,  X  =  0.4  , 

CL 

L  a 

=  4.0  , 

CAS  =  °-01  ■ 

CD0  " 

0.01  , 

s 

=  0.0625  , 

CD  =  CDo  +  sC l  =  0.0125  . 

Aircraft 

_3: 

CL  - 

0.8  , 

u> 

=  8.0  , 

k  =  0.4923  , 

v  as 

1  ,  X  =  0.4  , 

CL 

La 

=  4.0  , 

CAS  "  °-01  * 

CD0  = 

0.01  , 

s 

=  0.0625  , 

0^  =  Cj)q  +  sC^  =  0.05  . 

Aircraft  la: 

As  for  aircraft  1  but  with  s  increased  to  0.08,  which  gives  =  0.0612  . 

The  loci  of  the  coefficients  g,  h  are  traced  for  these  aircraft  when 
(a)  k  is  assumed  to  vary,  and  (b)  when  to  is  assumed  to  vary,  all  other  para¬ 
meters  being  held  constant  in  each  case. 
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Spitfire  V  (aircraft  of  R&M  2027)  : 

This  aircraft  was  used  in  the  numerical  examples  of  Ref  4  and  so  we  can 
quote  the  coefficients  of  the  stability  quartic  as  functions  of  m  and  k 
given  therein, 

B  =  6.010 

C  =  5.738  +  w 

D  =  0.1189  +  0.03039u  -  0.01763k 

E  =  0.002357w  -  0.07365k  . 

Static  and  manoeuvre  margins: 

K  =  0.003597m  -  0.1124k 

n 

H  =  0.01914  +  0.003432m  . 

m 

The  coefficient  of  m  in  has  been  adjusted  slightly  as  compared  with  the 

value  given  in  R&M  2027  in  order  to  render  <*  E  (as  it  should  be) . 

The  variation  of  k  with  fixed  values  of  m  and  variation  of  m  with 

fixed  values  of  k  are  considered.  An  analysis  in  terms  of  K  and  H  is 

J  n  m 

also  considered. 

Aircraft  of  R&M  2078  ^ : 

In  Ref  7,  Neumark  analysed  the  stability  of  an  aircraft  at  two  lift 
coefficients  in  terms  of  its  dependence  on  the  parameter  m  by  means  of  a  one- 
parameter  stability  diagram.  The  same  data  .re  used  as  a  further  example  of 
the  use  of  (g,h)  and  (G,H)  plots. 

In  this  case. 


CL  -  1,0  ’  CDo 

0.02  , 

CAS 

=  0.01  , 

s 

=  0.06  , 

3C 

L  .  c 

-t—  =  4.5  ,  v  = 

da 

3  , 

X 

=  1  , 

K 

=  0  . 

A. 4  Effect  of  varying  k  on 

the  stability  of 

aircraft 

1,  2 

and  3 

The  characteristics  of  the  three  aircraft  designated  aircraft  1 ,  2  and  3 
in  Fig  5  were  chosen  such  that  the  locus  of  g,  h  for  varying  k  corresponds 
to  a  different  one  of  the  three  cases  illustrated  diagramatically  in  Fig  1. 

Datum  values  of  h,  g  were  obtained  by  solution  of  the  stability  quartic 
(square  symbols)  and  also  by  use  of  the  relationship  of  equation  (A-4)  (circle 
symbols).  For  all  three  cases  the  locus  of  g,  h  is  close  to  the  straight  line 


x  lags SSSiSiiEBBg-  -■ 


4 
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approximation  obtained  by  treating  both  g  and  h  as  small  (see  section  3). 

In  the  case  of  aircraft  3  some  departure  from  near  linearity  is  evident  as  k 
is  decreased  and  the  region  of  unstable  oscillation  is  entered. 

In  the  case  of  aircraft  1  increase  of  <  takes  the  aircraft  out  of  the 

stable  region  and  just  into  the  region  of  unstable  oscillation  before  entering 

2 

the  region  where  the  factor  (X  +  gX  +  h)  corresponds  to  one  subsidence  and  one 
divergence.  If  the  value  of  30^/30^  is  increased  (aircraft  la)  this  trend  is 
more  pronounced  as  shown  in  Fig  6.  At  g  =  0  , 


K 


-  Ch  +  Rw 


T 


and  by  the  second  of  equations  (A-2) 


(A— 15) 


E  =  Ru)  -  Tx 


Ch  -  h2 


Thus  E  >  0  and  each  of  the  aircraft  1,  la  and  3  is  statically  stable  even  when 
k  is  such  that  it  enters  the  unstable  region. 

The  linear  approximation  to  the  locus  of  g,  h  for  varying  k  is 
obtained  by  ignoring  all  terms  on  the  right-hand  side  of  equations  (A-2),  that 
is,  equation  (A-4)  reduces  to, 

S  _  (Pv  +  Rx  +  Qui)  ~  (m  +  P  +  Nv  +  Qx)g 
T  Rw-(u)  +  P  +  Nv  +  Qx)h 


(Pv  +  Rx  +  Quo)  -  Cg 
Rw  -  Ch 


(A- 16) 


A. 5  Locus  of  g,  h  for  varying  m  ,  all  other  parameters  constant  for 
aircraft  1 ,  2  and  3 

We  now  examine  the  effect  on  the  values  of  g  and  h  of  varying  the 
parameter  w  for  aircraft  1 ,  2  and  3.  For  the  variation  of  w  in  the  neigh¬ 
bourhood  of  the  basic  aircraft  the  resulting  curves  correspond  to  three  of  the 
various  possibilities  illustrated  in  Fig  2.  At  g  =  0  , 


Bh 


(Pv  +  Rx  ~  Sk) 


h2  -  (P  +  Nv  +  Qx)h  -  Tic 
(h  -  R) 


(0 


or 


(A— 1 7) 
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and 


E  =  Ch  -  h2  =  (u>  +  P  +  Nv  +  QX)h  -  h2 

h[Rjh  -  (P  +  Nv  +  QX)f  -  Tk] 

_____  • 


(A- 18) 


The  curves  for  aircraft  1  and  aircraft  2  (see  Fig  7)  intercept  the  positive  h 
axis  at  h  =  0.15265  and  0.25064,  and  at  0.58798.  When  the  values  of  w  and 
E  are  determined  for  these  points  it  is  found  that  for  7.73  >  io  >  2.04  air¬ 
craft  1  enters  the  unstable  oscillation  region  even  though  E  (or  static  margin) 
is  positive.  Aircraft  2  enters  the  same  region  at  to  <  9.65  .  Aircraft  3  does 
not  exhibit  instability  for  the  range  of  values  covered  by  Fig  7,  which  corres¬ 
ponds  to  values  of  to  from  2.23  up  to  +°°  .  For  some  negative  value  of  to  it 

2 

is  possible  that  the  locus  of  g,  h  intercepts  the  boundary  g  =  4h  , 
of  Fig  1  or  8. 

The  exact  relationship  of  h  and  g  for  varying  to  (equation  (A-5))  can 
be  written  as  a  quadratic  in  h  (this  is  the  form  convenient  for  computation). 

If  we  write  the  quadratic  as, 

B2h2  +  B(h  +  BQ  =  0  (A-l 9) 

then 

B2  =  B  -  Q  -  g  , 

B,  =  |Q(P  +  Nv  +  QX)  -  (Pv  +  RX  -  Sk)  +  BR }  -  (BQ  -  2R)g  +  Qg2  , 

r\ 

Bq  =  |R(Pv  +  RX  -  Sk)  +  TQk|  -  {R(P  +  Nv  +  QX)  +  Txfg  +  Rg  (B  -  g)  . 


As  with  the  corresponding  equation  of  the  main  text  (equation  11)  it  is  possible 
to  approximate  to  this. 

Approximation  1 : 

Linearize  the  expressions  on  the  right-hand  side  of  equations  (A-3)  and 
the  following  results, 


g  -  Q  _  (Pv  +  Rx  -  S<)  -  Bh  -  (P  +  Nv  +  Qx)g 
h  -  R  -  Tk  -  (P  +  Nv  +  Qx)h 


(A-20) 
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which  yields  the  following  expressions  for  ,  Bj  and  Bq  , 


B  =  Q(P  +  Nv  +  Qx)  -  (Pv  +  Rx  “  Sk)  +  BR  , 

Bq  =  |R(Pv  +  Rx  -  Sk)  +  Tk}  -  |r(P  +  Nv  +  Qx)  +  T<(g  . 

Approximation  2 : 

Another  type  of  approximation  is  obtained  if  we  consider  both  g  and  h 

2  3 

in  equation  (A-19)  to  be  small.  In  the  simpler  of  these  terms  in  gh,  g  ,  g  , 

2  2 

g  h  and  gh  are  all  neglected,  then  the  expressions  for  B^,  Bj  and  Bq 
become, 

B2  =  B  -  Q  , 

Bj  =  Q(P  +  Nv  +  Qx)  -  (Pv  +  Rx  -  Sk)  +  BR  , 

Bq  =  |R(Pv  +  Rx  -  Sk)  +  Tk(  -  {R(P  +  Nv  +  Qx)  +  TK^g  . 

It  may  be  noted  that  the  coefficient  B2  differs  from  that  of  approximation  1 
and  that  approximation  2  is  somewhat  more  accurate.  A  further  approximation 
can  be  obtained  by  neglect  of  the  third  order  terms  only. 

Approximation  3: 

Yet  another  type  of  approximation  is  obtained  by  considering  g  only  as 

.  2  3 

small.  If  in  the  expressions  for  B0,  B^  and  Bq  terms  m  g  and  g  are 
neglected,  we  obtain, 

B2  =  B  -  Q  -  g  , 

Bj  =  |Q(P  +  Nv  +  Qx)  -  (Pv  +  Rx  -  Sk)  +  BR}  -  (BQ  -  2R)g  , 

PQ  =  Jr'Pv  +  Rx  -  Sk)  +  Tk|  -  {R(P  +  Nv  +  Qx)  +  TK^g  . 

In  this  fern  B2  is  exact  and  B|  is  an  improved  approximation.  If  only  terms 

in  g^  are  neglected  then  B|  takes  on  its  exact  form  and  Bq  only  differs 

3 

from  its  exact  value  by  -Rg 
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In  Fig  7,  approximation  1  is  used  to  obtain  values  of  h  and  g  for  com¬ 
parison  with  the  exact. 

A. 6  Effect  of  varying  k  on  the  stability  of  the  Spitfire  V 


In  Ref  4  the  coefficients  of  the  longitudinal  stability  quartic  for  the 

Spitfire  V  are  quoted  as  functions  of  gd  and  k  .  We  now  consider  the  locus  of 

g,  h  for  varying  k  at  four  levels  of  u>,  C  or  H  Thus  in  Fig  8, 

m 

id  =  -7.035,  =  -0.005.  With  this  negative  manoeuvre  margin  the  aircraft  is 

statically  stable  when  k  ^  -0.224  (K^  »  0,  see  Fig  8).  However  it  falls  within 

the  unstable  oscillation  region  for  values  of  k  <  -0.66  (or  K  >  0.05 

n 

approximately).  Between  K  =  0  and  0.05  the  motion  consists  of  two 

divergences.  In  fact  the  phugoid  motion  is  only  stable  when  the  static  margin 
is  more  negative  than  about  -0.007. 


Approximate  values  of  the  static  margin  are  marked  alongside  the 

curve.  As  remarked  elsewhere  variation  of  k  for  fixed  w  is  equivalent  to 

the  variation  of  K  for  fixed  H 
n  m 

The  next  figure  (Fig  9)  shows  similar  trends  for  the  case  of  the  aircraft 
with  zero  manoeuvre  margin.  With  manoeuvre  margins  of  0.01  and  0.05  the  locus 
over  the  range  of  of  -0.2  to  +0.2  does  not  depart  much  from  the  linear 

approximation,  see  Fig  10.  However  for  H^  =  0.01  and  Kn  48  -0.2  the  approxi¬ 
mate  values  of  g  and  h  are  appreciably  in  error. 


I 

1 


( 


Thus  far  we  have  been  concerned,  in  the  main,  with  the  locus  of  g,  h  in 

the  neighbourhood  of  some  chosen  datum.  To  each  point  of  the  curve  there 

corresponds  a  value  of  k  (and  hence  Kfi)  as  given  by  equations  (A-6)  and  (A-7) 

and  the  values  of  along  the  curves  of  Figs  8,  9  and  10  are  determined 

approximately  by  interpolation  of  the  k  values.  When  the  point  g,  h  (or  the 

corresponding  point  G,  H)  represents  an  oscillation  no  difficulty  arises  in 

the  calculation  of  k  .  However  when  all  roots  are  real  ambiguity  can  arise  as 

discussed  in  section  6.  To  illustrate  the  point  we  consider  that  part  of  the 

2 

locus  in  Fig  8  which  lies  below  the  curve  g  =  4h  .  This  is  reproduced  in 
Fig  1 1 .  Since  the  points  between  A  and  B  on  the  locus  of  g,  h  and  the 
points  between  the  corresponding  points  A'  and  B'  of  the  G,  H  locus 
represent  real  roots  of  the  quartic,  there  are  six  combinations  of  these  roots. 
These  appear  as  sets  of  three  points  on  each  curve  as  indicated  in  Fig  11. 

Apart  from  the  fact  that  the  values  of  k  must  be  in  sequence  as  the  chosen 
point  moves  from  A  to  B  ,  there  seems  to  be  no  other  way  of  discriminating 
between  the  alternatives  unless  a  one-parameter  stability  diagram  is  prepared 
(of  Figs  13  to  15). 


o 

O' 

Ln 
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It  is  of  interest  to  examine  the  trends  for  large  positive  and  large 
negative  values  of  the  parameter  k  .  For  numerically  large  values  of  g  we 
may  approximate  to  equation  (A-4)  as  follows, 


that  is 


2  2  3 

kh  -  kg  h  +  gJ  =  0 


(A-21) 


(A-22) 


The  solution  of  equation  (A-22)  is 

h  =  g2  or  &  .  (A-23) 

These  limiting  forms  of  the  curves  are  sketched  in  Fig  12. 

To  relate  the  value  of  g  to  the  trend  in  k  we  note  that  for  large  g  , 

gh  -*•  HSk  +  g3)  .  (A-24) 

When  equations  (A-23)  and  (A-24)  are  combined  we  have, 

g3  -*■  ±  Sk  .  (A-25) 

The  two  solutions  in  equations  (A-23)  and  (A-25)  refer  to  the  g,  h  locus  and 
the  G,  H  locus  respectively,  if  these  are  identified  with  the  'phugoid'  factor 
and  the  'short  period'  factor.  This  is  arbitrary  and  to  emphasise  this  the 
curves  of  Fig  12  are  marked  with  a  dual  interpretation. 

A. 7  Locus  of  g,  h  as  oo  is  varied  (Spitfire  V) 

If,  in  the  data  for  the  Spitfire  V  as  listed  earlier,  we  fix  k  at  -0.085 
the  coefficients  of  the  stability  quartic  become  functions  of  w  as  do  the  two 
margins  and  .  To  trace  the  locus  of  g,  h  as  u>  is  varied  equation 

(A-5) ,  or  its  equivalent  equation  (A-19),  is  used  to  evaluate  h  for  a  range  of 
values  of  g  .  The  locus  is  a  closed  curve  approximately  elliptical  in  shape, 
see  Appendix  B.  Along  the  arc  ABC  of  this  curve  the  values  of  u>  as  given  by 
equations  (A-8)  or  (A-9)  are  not  necessarily  unique,  as  it  turns  out  that  all 
four  roots  of  the  quartic  are  real  in  this  range.  Accordingly  there  are  three 
values  of  g  and  h  corresponding  to  a  given  w  in  the  range  affected  as 
illustrated  by  the  points  for  w  =  -3  and  -4  in  Fig  13a&b. 
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A  portion  of  the  locus  is  shown  enlarged  in  Fig  13b,  so  as  to  enable 
approximate  values  of  u>  to  be  marked  along  the  curve.  The  point  to  -  ±<®  is 
given  by  g  =  Q  =  0.03039  and  h  =  R  =  0.002357  .  Also  plotted  in  Fig  13a  is 
the  locus  of  g,  h  when  the  value  of  <  is  changed  to  -0.127.  It  also  passes 
through  the  same  point  for  to  =  ±°°  ,  as  would  the  whole  family  of  curves  for 
different  k  .  For  points  between  C  and  A  the  factor  X  +  Xg  +  h 
represents  a  stable  oscillation,  see  Fig  13. 

In  the  next  figure  (Fig  14)  the  locus  of  G,  H  (the  'short  period’ 

factor  coefficients)  is  plotted.  On  it  are  indicated  the  points  corresponding 

2 

to  A  and  C  of  the  curve  in  Fig  13  and  the  boundary  G  =  4H  (c/  Fig  1). 

The  points  for  u>  =  -3  and  -4  would  number  three  in  all,  corresponding  to 
different  combinations  in  pairs  of  the  four  real  roots,  of  which  only  one  each 
is  shown. 

The  minimum  value  of  G  corresponds  to  B  -  g  ,  where  g  is  the 
maximum  value  of  g  .  At  the  value  of  G  given  by  (B  -  Q),  u>  ±°°  and  the 
locus  becomes  asymptotic  to  the  line  G  =  B  -  Q  .  On  the  basis  of  the  results 
obtained  a  one-parameter  stability  diagram  may  be  prepared.  This  entails 
solution  of  two  quadratics  and  as  explained  previously  no  difficulty  arises 
when  either  quadratic  represents  an  oscillation.  However  between  u>  =  -4.914 
and  -2.60  the  quartic  has  four  real  roots  and  it  is  necessary  to  proceed  in  the 
following  manner  to  relate  these  to  a  certain  value  of  to  .  For  a  point  of  the 
g,  h  locus  (along  ABC)  determine  the  value  of  to  and  the  corresponding  value 
of  G  =  B-g  .  With  this  value  of  G  we  solve  for  H  and  accept  the  value  of 
H  which  yields  the  same  value  of  w  as  given  by  the  point  g,  h  .  Thus  the 
curves  of  Figs  13  and  14  contain  the  information  of  the  many  curves  of  Fig  15. 

As  to  -*■  I00  the  quartic  factorizes  approximately  as  {x^  +  (B  -  Q)X  +  to}  x 

2 

|x  +  QX  +  r|,  a  result  which  follows  easily  by  inspection  or  from  the  equations 
(A-5) ,  (A-8)  and  (A-9). 

A. 8  Locus  of  g,  h  for  varying  to  ,  all  other  parameters  constant,  for  the 

aircraft  of  R&M  2078 

In  Ref  1 1  Neumark  tabulates  the  roots  of  the  longitudinal  stability 
quartic  for  an  aircraft  of  specified  characteristics  over  a  wide  range  of  u>  , 

-2  to  The  results  are  also  displayed  in  the  form  of  a  one-parameter 

stability  diagram.  It  is,  therefore,  interesting  to  examine  the  locus  of  g,  h 
and  G,  H  for  this  aircraft  as  given  by  the  present  method. 
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In  Fig  16  the  coefficients  g,  h  as  derived  from  Neumark's  results  and 

further  values  as  obtained  by  application  of  equation  (A-19)  are  plotted. 

Neumark's  results  identify  the  point  A  as  corresponding  to  u>  *  0.2369  in 

agreement  with  a  value  of  w  =  0.2367  from  equation  (A-9) .  Since  he  calculated 

the  roots  for  the  two  negative  values  of  u>  =  -1.0  and  -2.0  on  the  arc  ABA'  it 

seems  that  he  might  have  attempted  to  locate  A'  ,  the  other  point  of  inter- 

2 

section  with  the  curve  g  =  4h  ,  but  failed  to  hit  on  the  appropriate  value  of 
w  .  Points  on  the  locus  above  A'  correspond  to  the  third  type  oi  oscillation 
arising  from  a  recombination  of  two  real  roots,  of  Fig  15.  The  locus  is  to  a 
close  approximation  an  ellipse  as  indicated  in  Fig  16. 

To  obtain  the  corresponding  G,  H  locus  equation  (A-19)  is  used  to 
calculate  H  for  a  range  of  values  of  G  in  the  neighbourhood  of  G  =  B-Q  , 
the  asymptotic  value  of  G  as  to  -*■  ±°°  .  When  the  coordinates  so  obtained  and 
those  derived  from  Neumark's  results  are  plotted  the  curve  of  Fig  17  is  obtained. 

A. 9  Analysis  of  the  stability  of  the  Spitfire  V  aircraft  in  terms  of  its 
static  and  manoeuvre  margin  values 

Since  the  static  margin  is  linearly  related  to  the  parameters  to  and  k 
and  the  manoeuvre  margin  to  the  parameter  to  ,  we  may  write  the  coefficients  of 
the  quartic  in  terms  of  the  margins  and  ,  as  follows, 

B  =  6.010  , 

C  =  291 . 37529H  +  0.16108  , 

m 

D  =  8.69050H  +  0.I5685K  -  0.04744  , 

m  n 

E  =  0.65525K  . 

n 

These  enable  the  stability  of  the  aircraft  to  be  examined  using  equations  (A— 1 3 ) 

and  (A-14)  when  K  is  varied,  for  fixed  H  and  H  is  varied,  for  fixed 
n  mm 

,  respectively. 

A. 9.1  Locus  of  g,  h  for  varying  K  ,  fixed  H 
_ n _ m 

As  already  remarked  the  variation  of  K  with  a  fixed  value  of  H  , 

n  m 

the  manoeuvre  margin,  is  equivalent  to  variation  of  k  with  to  held  fixed. 
Accordingly  the  curves  of  Figs  8  to  10  are  marked  with  the  values  of  and 

the  fixed  value  of  Hm  indicated  on  each  diagram  and  are  capable  of  the  dual 
interpretation. 


*■ 

i 


32 


Appendix  A 


A. 9. 2  Locus  of  g,  h  for  varying  H  ,  fixed  K 
_ m _ n 

To  fix  the  value  of  ,  the  static  margin,  is  to  establish  a  linear 
relationship  between  to  and  k  .  Hence  as  (or  to)  is  varied  the  value  of 

<  changes  to  retain  the  constant  ,  static  margin.  With  =  0.05  the 
locus  of  g,  h  for  varying  as  given  by  equation  (A-14)  is  the  curve  shown 

in  Fig  18.  It  follows  from  equation  (A-14)  that  as  h  -*•  ±0,  g  -*■  0.0298  and 
using  the  associated  equations  for  Hm  it  can  be  seen  that  H  -*•  ±»  . 

By  use  of  the  same  equation  it  is  possible  to  construct  the  locus  of 
G,  H  ,  the  coefficients  of  the  other  quadratic  factor  of  the  stability  quartic. 
This  is  illustrated  in  Fig  19  where  it  can  be  seen  that  G  =  B- 0.0298  is  an 
asymptote.  This  corresponds  to  the  point  0.0298,  0  of  Fig  18.  The  result  can 
be  readily  derived  by  approximate  factorization  of  the  quartic  as  -*■  °°  . 

The  quartic  becomes  approximately 


X4  +  BX3  +  C  H  X2  +  D  H  X  +  E  K  =  0  , 

mm  mm  n  n 


which  factorizes  into 


* 

« 

•  1 

_  /  D  \ 

_  D  E  K 

X  +  C  H 

►  - 

x2  +  _m  X  +  rnun 

\  c 

m  m 

C  C  H 

\  mJ 

m  mm 

.  . 

>  =  0 


That  is. 


B  -  7T-  ,  H  ->•  C  H  -*•  ±  «> 

C  mm 

m 


and 


D 

m 

C  ’ 
m 


h  -*■  ±  0  . 


The  curves  of  Figs  18  and  19  cover  the  range  of  practical  interest,  but  the  fact 
that  curves  are  not  bounded  in  the  same  way  as  the  curves  for  varying  u  is  of 
interest  and  merits  further  consideration. 

Equation  (A-14)  can  be  written  as 


+  Cjh  +  C 


0 


0 
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where  C2  «  B-g  - 

m 


g2  -  B  —  g  +  I?  +  (N  -  iCLc()v  +  QX ( 
m  m  m 


D 

'<T 

i  m 


For  large  values  of  g  the  coefficients  C2>  Cj  and  Cq  are  given  approxi¬ 
mately  by, 

C2  »  -8 
Cj  ~  cg(g  -  B) 


c0  -  Eg  , 


where  c  *  D^/C^  •  This  gives  a  quadratic  to  which  the  solution  is, 


h  98  c(g  -  B)  or 


c(g  -  B) 


Thus  in  the  limit  h  behaves  as  either  eg  or  -(E/cg)  (E  >  0,  c  >  0)  .  If 
we  identify  g,  h  with  a  point  on  the  branch  of  the  curve  of  Fig  18  that  refers 
to  the  'phugoid'  factor,  then  we  are  interested  in  the  behaviour  of  h  as 
g  -*■  .  It  follows  from  the  result  just  given  that  h  -*■  -®  or  +0  ,  and  that 

-*•  -°°  in  both  cases.  When  g  -*■  -<*>,  G  ->  +»  since  G  +  g  =  B  .  By  virtue  of 
the  solution  to  the  quadratic,  which  applies  equally  to  G  and  H  ,  we  have 
H  -*•  +«,  -0  ,  whilst  again  H  •  The  curves  of  Figs  18  and  19  show  that 

for  Hm  -*■  the  quartic  has  four  real  roots  and  so  coefficients  g,  h  and 
G,  H  appear  three  times  on  each  curve  for  H  -*■  -»  ,  corresponding  to  the  six 
different  combinations  of  the  roots  in  pairs.  It  is  possible  to  derive  the  same 
result  by  direct  approximate  factorization  of  the  quartic.  For  large  Hm  the 
quartic  is  approximately  , 

in 

i  x4  +  bx3  +  c  h  x2  +  d  h  x  +  e  = 

r-  m  m  m  m 


0 
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which  factorizes  into 


0 


where  C  H  =  -k^  , 
m  m 

quartic  as  H  -*•  -« 


k  large  and  positive, 
are 


Thus  the  four  real  roots  of  the 


(a)  A 

(b)  A 

(c)  A 


k  +  IUH-B 


2  C 


m 


m 


m 


(d)  A 


EC 

m 


approximately  . 


If  we  assume  k  positive, 

(a)  and  (b)  yield  the  point  G 


(b)  and  (c)  yield  the  point  G 


(b)  and  (d)  yield  the  point  G 


(a)  and  (c)  yield  the  point  g 


(a)  and  (d)  yield  the  point  g 


(c)  and  (d)  yield  the  point  g 


* 

m 


-+  —  no 


kD 


k  ;  H 


m 


m 


k  -*•  +  °°  ;  H 


EC 

m 

D  k 
m 


-+•  -  0 


-  k  -  oo  j  h 


D  k 
m 

C 

m 


—  oo 


-  k  -*■  -  oo  j  h 


EC 
_ m 

D  k 
m 


-*■  +  0 


m 

C  * 
m 


E 

7 


This  explains  the  trends  exhibited  by  the  curves  of  Fig  20. 
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The  difference  between  the  closed,  compact  nature  of  the  curves  of 
Figs  13  and  16  as  compared  with  that  of  Fig  20  underlines  the  essentially 
approximate  nature  of  the  manoeuvre  margin  as  a  system  parameter.  In  contrast, 
the  concise  quantity  u  is  a  natural  system  parameter. 


f 


1 
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SOME  REMARKS  ON  THE  GEOMETRY  OF  THE  LOCUS  OF  g,  h 

A  look  at  the  approximations  to  the  various  loci  to  a  second  order  in  g 
and  h  (the  resulting  curves  are  then  conics)  helps  to  establish  the  nature  of 
the  curves  generally. 

B.  1  Locus  of  g,  h  for  varying  <  (second-order  approximation) 

From  equation  (A-4)  we  have,  to  second-order  in  g  and  h 


(bH 

(B- 


gh  -  Bg  +  (w  +  P  +  Nv  +  Qx) g 


~  (to  +  P  +  Nv  +  Qx)  >  h  + 


(»  I  -  «)■ 


+  Pv  +  RX 


The  nature  of  this  conic  is  determined  by 


1  I  b2s2 

4  I  t2 


+  4  <  0 


Hence  the  curve  is  a  hyperbola  to  this  degree  of  approximation. 

B.2  Locus  of  g,  h  for  varying  oj  (linear  approximation  of  equation  (A-20)) 

The  equation  to  the  curve  according  to  approximation  I  (see  equation 
(A-20))  is 

Bh2  +  |Q(P  +  Nv  +  QX)  -  (Pv  +  RX  -  Sx)  +  BR}h 

+  |R(Pv  +  RX  -  Sx)  +  Txj  -  {R(P  +  Nv  +  QX)  +  Tx(g  =  0  . 

.  2 

Since  there  is  no  term  in  gh  or  g  this  is  a  parabola.  Furthermore  since 
the  equation  to  the  curve  can  be  written  in  the  form. 


B(h  -  hQ)2  -  |R(P  +  Nv  +  QX)  +  Tic }  (g  -  gQ)  , 


the  axis  of  the  parabola  is  parallel  to  the  h  *  0  axis. 
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B. 3  Locus  of  g,  h  for  varying  id  (quadratic  approximation  of  equation 
(A-20) ) 

The  equation  to  the  curve  in  this  case  is 


+ 


Blh 


+  Br 


0  , 


where 


B  -  Q  , 

Q(P  +  Nv  +  QX)  -  (Pv  +  RX  -  Sic)  +  BR  -  (BQ  -  2R)g  , 

2 

{R(Pv  +  RX  -  Sk)  +  Tk|  -  |R(P  +  Nv  +  Qx)  +  Tic}g  +  BRg  . 


This  is  the  next  higher  order  approximation  mentioned  in  Appendix  A,  section  A. 5 
(approximation  2) . 

2  2 

Here  the  coefficients  of  h  ,  gh  and  g  are  as  follows, 


h2  :  B  -  Q 
gh  :  2R  -  BQ 
g2  :  BR  , 


and  the  test  function  for  the  conic  is 


But, 


BR(B  -  Q)  -  ^  (2R  -  BQ)2 


R 


v  8Cl\ 
+  2  av ) 


Q 


+ 


2  av/e 


2 

If  the  change  with  speed  of  R  and  Q  is  negligible  or  absent  then  R  =  , 

Q  *  Cn  and  the  test  function  for  the  conic  becomes, 
e 


4 


which  is  usually  >  0  ,  so  that  the  curve  is  an  ellipse. 
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LIST  OF  SYMBOLS 


Main  text 


Z 

C 


D 


C 


L 


D 


g 

G 

h 

H 


iy 

iy 


L 


m 

e 


in 


Q 

R 

Q  with  suffix 


u 

U  with  suffix 
V 


w 

W  with  suffix 


mean  chord  of  wing 

2 

drag  coefficient  =  D/jpV  S 
lift  coefficient  =  L/JpV^S 
drag 

coefficient  of  A  in  quadratic  factor 
coefficient  of  A  in  quadratic  factor 
constant  term  in  quadratic  factor 
constant  term  in  quadratic  factor 
manoeuvre  margin 

inertia  parameter  =  ly/n^c^ 
moment  of  inertia  about  y-axis 
coefficients  of  the  stability  quartic 

static  margin 

lift 

aircraft  mass  in  equilibrium  flight 


concise  moment  quantities,  see  R&M  3562  or  ESDU  data  items 


2 

kinetic  pressure  =  $pV 
Routh's  discriminant 

polynomial  involved  in  coefficients  of  stability  quartic 

component  of  aircraft  velocity  along  x-axis 

polynomial  involved  in  coefficients  of  stability  quartic 

resultant  velocity  of  aircraft 

component  of  aircraft  velocity  along  z-axis 

polynomial  involved  in  coefficients  of  stability  quartic 
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LIST  OF  SYMBOLS  (continued) 


concise  force  quantities,  see  R&M  3652  or  ESDU  data  items 


angle  of  attack 

root  of  the  stability  quartic 

relative  density  parameter  =  me/jpeSc 


air  density 

pendices  A  and  B  (of  R&M  2027  and  R&M  2078) 


B,C,D,E 


coefficients  of  stability  quartic 

propeller  contribution  to  xu  derivative  =  [ (-1 /pSV) (31/ 9V)]  , 

propeller  'drag'  coefficient 

rate  of  change  of  with  a  =  9C^/9a 

2 

inertia  parameter  =  Ig/me^ 
moment  of  inertia  about  y-axis 

tail  arm  (distance  from  eg  to  mean  ;-chord  point  of  the 
tailplane) 


+  *CD  +  CAS>< 


(CD  +  CAS>  (CLa  +  V  +  CL  + 


1(CD  +  CAS}  ~  *CL  tan  Yle 


$  k  -  m 


[jCL{CL  *  !V  TT  -  <CD  *  CAS>  ’'] 

i(^). 

*  c»  *  (cl  -  T?)  v} 


LIST  OF  SYMBOLS  (concluded) 


Note:  In  these  expressions  the  contribution  of  the  Z  force  due  to  rate  of 

pitch  has  been  neglected,  of  Table  2  of  R&M  2078.  In  the  form  given  they  apply 
to  a  propeller  driven  aircraft  at  speeds  for  which  compressibility  and  aero- 
elastic  effects  are  negligible.  In  the  general  case,  where  compressibility  and 
aeroelastic  effects  are  also  included,  the  C^g  term  needs  to  be  replaced  by 
CAS  +  *v(3Cd/3V)  *  cAs  is  zero  for  a  jet-propelled  aircraft. 

Y  angle  of  climb 

U  relative  density  parameter  =  m^/p^Si 

K 

concise  quantities  in  the  old  notation  (for  relation  to  concise 
'  ,  quantities  in  the  new  notation,  see  R&M  3562,  Part  5.  For  definition 

X  in  terms  of  old-notation  derivatives,  see  Refs  A,  5  or  6) 
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one  divergence 
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Fig  1  Diagram  illustrating  analysis  of  stability  using  locus  of  the  coefficients 
g,  h  for  varying  mu 

/  increasing  mw 
direction 


O  9  =  Wm2 

h  =  Wmi 


mw«+( 


Fig  2  Diagram  illustrating  analysis  of  stability  using  locus  of  the  coefficients 
g,  h  for  varying  m... 


Fig  5  Examples  of  g,  h  locus  for  varying  k  for  aircraft  1,  2  and  3  of  Appendix  A 


\ 


increasing  Kn 
K>r-JC)direction 


-  0*04 


basic  aircraft 


Fig  8  The  g,  h  locus  for  varying  Kn  or  k  (w  fixed  to  give  Hm  -  -0.006) 
(Spitfire  V.  Ref  6)  m 


-  00 1 


□  exact 


increasing  Kn  ^  basic  aircraft  X  approximate 
direction 

iU  rha  g,  h  locus  for  varying  Kn  or  k  (Hm  -  0.01  and  0.06) 

(Spitfira  V,  Raf  6) 


Fig  13a&b  The  g,  h  locus  for  varying  u>  and  fixed  x  (Spitfire  V)  showing 
multiple  points  which  occur  for  =  -3  and  -4 


4h 


8l 


o  Points  corresponding 
to  oj  values  for  which 

A2  +  gA  +  h  has  equal 
roots 


Asymptote 


oo  =-io 


CO  ■  -14 

CO 

* 

-CO 


Fig  14  The  locus  of  G,  H  for  varying  to  (x  =  —0.085)  associated  with 
g,  h  locus  of  Fig  13  (Spitfire  V) 
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Fig  15  A  one-parameter  stability  diagram  constructed  from  the  results  of 
Figs  13  and  14 


Fig  16  Locus  of  g,  h  for  varying  u  (aircraft  of  Ref  7)  with  comparison  of  results 
obtained  by  present  method  and  solution  of  quartic 
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Fig  17 


i 


Fig  18 


Fig  18  Locus  of  g,  h  for  varying  manoeuvre  margin  and  fixed  static  margin 
(Spitfire  V) 


Fig  19  Locus  of  G,  H  associated  with  that  of  Fig  18  for  varying  manoeuvre 
margin  and  fixed  static  margin  (Spitfire  V) 
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CORRIGENDA 

Page  4,  section  2,  line  1:  R&M  3562  is  Ref  10. 

Pages  9  and  10,  section  3.1:  The  equations  in  this  section  should  contain 
the  factor  not  (K^-k^)  as  printed. 

Page  13,  section  5:  The  table  of  aircraft  characteristics  given  at  the 
foot  of  page  13  and  top  of  page  14  is  incorrect.  It  should  read  as  follows: 
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Page  14,  penultimate  paragraph  of  section  5  is  incomplete.  It  should  be 
amended  to  read: 

"...  becoming  divergent,  when  condition  (i)  of  page  10  is  fulfilled.  If 
condition  (iii)  applies  increase  of  m^  achieves  the  same  result." 

Page  14,  last  sentence  of  section  5:  The  value  of  g  quoted  should  be 

0.06. 
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A  more  complete 
from  50  to  15  and 


version  of  Fig  4  is  attached.  This  corrects  the  value  of 
inserts  additional  values  along  the  curve. 
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